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ABSTRACT 



An analysis is done of the motion of an underwater cradle for submerged 
recovery of the DSRV Alvin. 

The cradle is modeled as a three degree of freedom system with four cor- 
ner springs and includes heave pitch and roll motions. The surface ship Lulu 
is taken as the driving force that couples energy into the cradle through the 
springs. Hydrodynamic damping values are found by a 64 section finite ele- 
ment analysis and are compared to corresponding linear approximations. 

Solution of the resulting equations was accomplished using a Runge-Kutta 
numerical integration method developed at M.I.T. 
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IPf will substantially decrease the transmit Ped ship motion. Additionally, 
the results indicate that the pitch and roll natural frequency is higher than 
the heave natural frequency and consequently is the limiting design considera- 
tion. 



The results of the analysis are discussed and several possible approaches 
are considered for producing a soft spring constant while at the same time 
supporting the high static load of the cradle weight. 
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NOiM}'I!CLATURE 



Sy:nbol 


Interpretation 


X 


Forward direction of Lulu 


y 


Downward motion of Lulu 


Z 


Athwart ships of Lulu 


^0 


Maxiraum heave amplitude of Lulu 




Maximum pitch angle of Lulu 


^0 


Maximum roll angle 


e 


Cradle pitch angle 




Cradle roll angle 




Heave and pitch frequency of Lulu 


^2 


Roll frequency of Lulu 


p 


of bccLwciucj. 




Coefficient of viscous damping 


Mo 


Mass of cradle 


Ma 


Added mass of cradle 


M<p 


Total mass of cradle 


M 

O 


Actual pitch inertia of cradle 


^al 


Added pitch inertia of cradle 




Total pitch inertia 


IQ2 


Actual roll inertia of cradle 


Ia2 


Added roll inertia of cradle 


^2 


Total roll inertia of cradle 


L 


Length of cradle 


V7 


Width of cradle 
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Tv. .^'"1 


Tnterpreta' ion 


Percent 


Percent deviation in springs 


:<i - :<!+ 


Four cradle support springs 


a 


Fore and aft sway 


B 


Transverse sway 


“l 


Coefficient of real part of complex pole 




Coefficient of imaginary part of complex pole 




Bulk modulus 




Flow to load 


Pm 

m 


Load pressure 


R 


Proportional control constant 


S 


Integral control constant 


£ 


Area of piston 


V 


Volume of piston 


Ke 


Cylinder compliance 


Ket 


Total compliance 


Kox 


Flow sensitivity of valve gain 


Kop 


Negative slope of Pjjj versus curve 


Fi 


Valve forces 


t 


time 


c(t) 


Set point input to P I control 


D 


Damping factor 


^11-^hk 


The sixteen elem.ents m.aking up quadrant one of the 
cradle's projected area 


6 


Spring stretch 
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CHAPTCR I. 

INTRODUCTIOIl 

The handling of heavy objects through the air-sea interface is not a 
nev.' or unique problem and has been studied many times in the past. The 
deep s'jbmergence research vehicle (DSRV) Alvin owned by the U.S. Navy and 
operated by 'Noods Hole Oceanographic Institute (VJHOI) is most likely the 
DSRV with the world's best record for successful dives completed. The 
submersible is handled by its support ship, a catamaran named Lulu and is 
transported to and from the open ocean dive site between Lulu's twin hulls. 

The launching and retrieving operation is accomplished with a great 
deal of care, Alvin being initially positioned on her 30-ton capacity cradle 
just aft of amidships. The cradle, (Figure 1), is supported at its four 
comers by chain hoists which are two-blocked such that the cradle is at 
deck level during transit. This allows easy maintenance, repairs and man- 
ning of Alvin . 

During launching, the catamaran lies with its stem to the seas, the 
cradle is then lowered and the floating submersible clears Lulu with the aid 
of Alvin's screw and nimerous snubber lines. Recovery is the reverse of the 
operation, Alvin surfacing well away from Lulu, closing Lulu, and finally 
driving between Lulu's hulls and being manually positioned over the sub- 
merged cradle. Last, when the wave timing is correct, the cradle- 
submersible combination is raised clear of the water. 

This method of launching and retrieval has allowed reasonably routine 
operation^ to occur in sea states up to three , with occasional sea state 
four retrieval due to emergency conditions. 



10 




Figure 1 
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V/ilhout goin.f. late a discuasicn of sea state probabilities, it is v/ell 
kncv;n that along the north Atlantic off Cape Cod the sea state vjill frequent- 
ly go from state three to state four due to a prevailing southv/est wind which 

increases from ten knots in the morning to about 20 in the afternoon. Sea 
state four makes for a very difficult recovery should it occur during the 
time of a dive. The problem here really is that once the Alvin is in posi- 
tion to be hoisted aboard, the relative motion betv;een Alvin, the cradle, and 

Lulu could cause very high and damaging shock loading on Alvin, the hoist, 

or the cradle. This loading could possibly result in severe damage to Alvip , 
or even another loss. 

Clearly then, a smooth docking to a stationary cradle followed by a 
rapid hoisting through the air-sea interface would be highly desirable and 
won! a vT-ntuailv eliminate shock loading and consequent damage. 

Fortunately, the Alvinas cradle consists of four 120* lengths of chain 
and four synchronized wildcats. If the cradle were lowered to approximately 
100* below the surface, the surface wave effect would be virtually inefec- 
tive on Alvin and it would be possible to dock or depart from the cradle 
without the fear of poor surface conditions. The problem which appears to 
be the limiting factor is the motion of the cradle itself. This motion is 
nearly the same as that of Lulu due to the Inelastic chains which connect 
the two. Some work has been done to show the transmitted surface motion 
that the cradle would feel. It has been generally concluded that if members 
more elastic than chains were used, perhaps the cradle motion could be un- 
coupled from the Lulu and the Alvin could dock with a virtually motionless 
cradle at 100* and then be hauled to the surface rapidly and safely through 
the air-sea interface. 
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Vao corrplote -^.na lysis of the Lulu-cradle interaction is the primary 



ncern in this thesis. 
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CHAPTEK II. 

Ti-iEOPY a::d development of eouatio;;s of motion 



To this tine the ship cradle interaction motion has been studied using 
a simple model shov;n in Figure 2. While this model applies very well to a 
svstem v;ith a single support it does not predict the motion of a four point 
suspended cradle with ship pitch, roll, and heave. 




M 



A 

X^=Yosir.wt 



I 

I 

I y(t}=yQsino)t 




D 
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Figure 2 . 

Looking at Figure 3 it is apparent that the motions of a cradle be- 
neath a ship are extremely complex. The cradle can exhibit vertical motion 
(y), pitch (9), roll (d), fore and aft sway (a), transverse sway (B), twist- 
ing about the y axis, as well as horizontal motion in the x direction. In 
efrect, all six degrees of freedom can be achieved with the cradle. For- 
tunately, to somewhat simplify the problem, fore and aft sway (a). 
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Figure 3. 
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transversc sway (g), ano tv;ist v;ere c:<porimental.1y found to be insignifi- 
cant . 

Consequently, the problem to be solved is a three degree of freedom 
one v:ith roll (cj)), pitch (0), and heave (y) , and where all three modes of 
motion can transfer energy back and forth thus contributing to the com- 
plexity of the problem. 

Before beginning the modeling the initial conditions must be defined. 
To begin with, the cradle is assumed to be of length 1 = 17.5 feet, width 
w = 16.3 feet, with a v/eight of approximately 508 slugs. Since the 
reference material has indicated that a possible mode of operation would 
be with a flat plate affixed to the bottom of the cradle, this is the mode 
chosen for our model. Consequently the cradle is assumed to be a rectan- 
gular flat plate with mass equivalent to that of the cradle alone. In the 
final analysis it will be shov/n that the added mass terms associated with 
the cradle so overpower the actual mass that the exact configuration of 
the cradle structure is inconsequential. 

It was first necessary to pick an adequate co-ordinate system which 
would accurately describe the system and make it possible to completely 
position the cradle knowing only three co-ordinates — heave (y), pitch (0), 
and roll . A reference system (Figure 4), fixed above the cradle was 
chosen with X positive for^ward, Y positive down and Z positive to port. 
Using the conventional right hand rule , positive 6 and positive ((> are as 
shown in Figure 4. As is conventional in dynamic analysis the cradle is 
given an instantaneous displacement positive y, positive 6, and positive ((> . 

Since the equations of the system which are being dealt with are 
quite complex with many terms, it will be convenient to set up the 
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Figure H. 



Cradle Equilibrium 
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rosultant equation in standard state variable form consisting of six first 
order differential equations and an output expression. The most general' 
form of such an equation is 

Xj^(t) Aj_ 2 Xq(t) ^>11^12 • • • ’^In U^Ct) 

X 2 (t) A21 A22 X 2 (t) 

+ 

^nl ^2 ^nn l^nl ^nr 

X(t) = AX(t) + Bu(t) 

where A is the N x N system matrix with B an N x r control matrix, 
btrictly speaking, tnis above representation applies only to a m •" oraer 
linear equations , however, since there will be a multitude of terms, both 
linear and nonlinear, it will be a good bookkeeping aid to arrange terms 
as such. . 

The basic equation of motion for a spring, mass, and damper system 
is used here. Development of the equations of motion is begun with the 
consideration of this free body diagram. 



X„<t) 

or 
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where : Z forces = ma 

or: 



t 



WEIGHT 



Figure 



" ^drag “ ^ spring + force (exciting) 
and is the case of three degrees of freedom 

“ *” ^ ~ -t” c 'nr'i n cr t -tAynitlnCX VY./ 



^drag ^spring ^exciting 

^2^ = - ^drag “ ^spring ^exciting ^^0^* 

Each term will be discussed independently. 



The Spring Force 

The four springs (chains) at the comers of the cradle are the major 
factor in decoupling the surface motion from the cradle. The spring force 
is proportional to its stretch and this stretch may well be different for 
each of the four springs. This is one of the mechanisms by which energy 
is transferred from one mode ^ such as heave, to pitch or roll. Ideally, 
if all four springs were equal and the cradle were given a pure pitch, roll, 
or heave motion, then the other modes v/ould not be excited and you would 
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havG sirnp-ly uncouDled systems, ilowever, considering the inexactnes 

of the Peal V/orld, the cross coupled damping, and the exterior forces, it 
v;ould be unrealistic to perform such a simple analysis. In actuality the 
ship pitches, rolls, and heaves about an equilibrium position as does the 
cradle. The stretch of the four springs is defined as the difference in 
position betv:een the position of the ship and the position of the cradle 
relative to their respective equilibrium for each motion mode. 

The respective stretches of the four springs K2 , K3 , and are 
as follows : 

= (y-YQsinuj^t) + (L/20-L/20Qsinu)j_t ) + (w/2c{)-w/2^Qsinm2t) 

.62 = (y-Y^sincoj^t ) + (L/20-L/20Qsino)2^t ) -- ^'w/2<{>-w/2^Qsincj2t ) 

63 = (y-rQSincjj^t) - (i,/2^-L/20^sini^}jt) - « w/2cp“W/2<PQSina32t; 

^4 = (y“YQsino3j_t) - (L/ 20 “L/ 20 ^sina)j_t ) + CW/2(J)“W/2<^>sinco2t ) 

This is assuming thatYQ,0^, and , are the ship’s pitch, roll, and 
heave. 

Carrying this a bit further, when the stretches are multiplied by 
their spring constants, we get the respective restoring spring forces and 
moments in general form. 

The equation for heave — 

my = _ - K2«2 - KgSg - 

The equation for pitch — 



lj_0 - - ” ^2^2 ^ ^4^4 
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The eouation for roll-- 



’* W V/ W ^ 

1^4, = + K 2 - 62 + K 3 li 53 - J 5 i - K 4 - 6 ^ 



or plugging in 6’s, 
heave equation-- 



iny ^2^2 **' ^3^3 ^14,^11 



Plugging in: 



^2 = 



K^C(yi 



- YQsinu^t) + (- 73 - 



K 2 C( 73 ^ - Yosinwj^t) + (^ 73 - 

L 

K3C(7i - Y^sinoj^t) - (j 73 - 
K4[(7jl - Y^sinuj^t) - ( j 73 - 



L „ . s /W 

2 ^2 ^5 

- eoSintUjLt) - (-75 

7 ^ j ys 

- e^sincoit) + (y 75 



w 



J $Q sinta2t)] 

j $Qsino)2t)] 

j $QSino)2t)J 

W , . NT 

y 4 >QSino) 2 t) J 



pitch equation-- 



IiYh = 



L L L W W 

+ K3 ^(7 jl - Y^sinw^^t) - (jY3 - 2 “ ^2 ^5 ~ J ^ 



+K, 



4 |t(yi - Yo^inwj^t) - (y 75 - ^ eQsincoj^t) + (^ 75 - j $Qsinco2t)] 



-Kjl ^(y^-YQsinwjLt) + 73- SpSimoj^t) + ( J 75- ^ $Qsincii2t)] 

L Jj jj 

"^2 ^(Yi-Yosinwit) + (-73- - "^7^5” y 



roll equation-- 



l 2 Y 6 = ■ Kl^^Yi-YoSincoj^t) + (y 73 - y egSinw^L''^) + (| 75 - y 4 >oSinw 2 ^n 



! 



I 

I 

( 
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^ ^ J ^5" f ">'oSinoj2t)] 

+K2^-{(yi-YoSinw^t) t (^ Y3- “ 9osinUj,t) - ij Y^- | '3>Qsimo2t)] 






Driving Force 

The force v/hich excites the cradle is transferred through the springs 
to the cradle. In this analysis it is assuir.ed that Lulu is heaving and 
pitching sinusoidal with the same frequency co]^, and amplitudes and 0 q , 
respectively. Additionally, it is assumed that Lulu is rolling sinusoidal 
ly v;ith amplitude and frequency o)2. 



Drag Forces 

Triis is the most difficult set of terms to analyse since due to the 
low Reynold's numbers involved the drag forces are primarily caused by 
hydrodynamic form drag rather than viscous drag and are nonlinear. In ad- 
dition, for geometric shapes with hard edges the drag coefficients are 
generally independent of Reynold's numbers. 

The general equation describing drag forces on such a shape is given 
as Drag = Cd V , and it is in this where we find the problem. If 

drag is to alvzays oppose motion as it does here , the equation above should 
actually be written Drag = Hq A^Cd|v|v, vzhich vzill preserve the correct 



sign. 






I 

I 
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Init is.lJ.y 8,n 6X5.Co c.'ti'Cciup'C o/z 'ons oGj.iibioxi Vc,5 'bi'icd., Tns.3iii\ici* c.3 3. 

differential element of drag is defined in quadrant (l) of Figure 6 as 

d(drag) = ^pCjjdxdz (y+x6cos4>cos0+Z(}>co3(iicos9 ) | y+x9cos4icos6+z(})cos(|)cos0 | . 

then the total drag on quadrant (l) would logically he 

DRAG = / / ^pC-p)dxdz(y+x9cos9cos<{i+z6cos0cos<}) ) |y+x0cos0cosi})+z$cos0cos4i I . 

0 0 



This approach could he pursued to an exact solution if an identity for 
(A+R+G) ! A+3+c| was obtainable. However, (A+B+C) |a+B+c 1 only equals 
(A+R+G) (A+3+C) when (A+B+C) is a positive number. The determination of 
■-he exact sign of (A+B+C) or (y+x9cos0cos(})+zit)cos9cos<}) ) can only be obt.^ined 
wr.''n a hr.cwn value of X, or Z is available. Consequently, if X and Z are 
tr.e i.ntegrand they cannot be assigned a specific value as they vary fro.m 






f 






I 
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0 to L/2 oP.'l V//2 r^jLipocti vely . Thus, as ina Lheir.atically satislying as c.s 
exact equation for drag v;ould be , it is alas, impossible to obtain. Due 
this, a more straight forward yet more tedious finite element method v;as 
used to approximate drag forces. 

The cradle was first divided up into four quadrants as shown below: 



to 




Figure 7. 



Using the sign convention previously described, the velocities of any point 
in the four quadrants are as follows : 

Vi = (y + X 0 cos0cos({) + z({^cos0cos({)) 

V 2 = (y + X 0 cos0cos({) - Z({)cos0cos4) ) 

V3 = (y - X 0 cos0cos4> ~ Z(|)cos0cos(J)) 

^4 ~ ^ ® COS0COS(() + Z(f)COS0COS({) ) . 

Each quadrant is then divided into sixteen areas and the drag is calculated 
for each one as shown for element below. 




'■h 





-?.h- 



Dra.^T olcaient aj 



22 



= h. PCj;) Area (&22^ ^^221^^221 



. W L 

Area a 22 = “ ' g 



V 22 = (y+ L0 cos0cos!})+ V/({icos0cos(}) ) 

16 ID 



2 



ail 


‘'^12 


ai3 


ail, 


^21 


^22 


^23 


^2k 


^31 


^32 


®-33 


B-2k 


aiil 


ai^2 


^h3 


^hk 



JJ2 



( Quadrant 1 ) 
Figure 8. 

W L * 3_ • 3 • I • 3 • 

Drag(a 22 ) = ^ '^^0cos0cosi})+ ^//ij)cos6cos<f> ) |y+ |j^L0cos0cos(j) 



3 • I 

+ ■^/f(|)COS0COS(f)] 



Hence, the velocity is assximed to he fixed over the entire element and thus, 
so IS zne sign oi une arag. laeaixy , as tiie xiumuei’ ui t?j.cmcxiot> uj. , ihc 

possibility of the total velocity changing sign in an element goes down un- 
til, in the limit, you have numerically integrated the whole surface to an 
exact result. In this analysis 64 elements are taken each time step and 
added together to produce the resultant drag. 

Pitch and roll drag are calculated in much the same way hut with the 
added complication that each drag element acts through a finite moment arm 
to produce the drag in terns of foot Ihs. One difference is the proper 
description of the velocity of the element which for our same element a.22 
would he 



Vapp = (ycos0cos(j) + ~7 L 0+ 

16 






The computer subroutines 'for calculating the drag forces are given in 



the Appendix. 
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In cases where the relative velocity betv;een a submerged body and its 
surrounding water does not remain constant, additional forces come into play 
as a result of the acceleration of the mass of water entrained by the body. 
This added mass depends on the physical dimensions of the body as well as 
its mode of motion. The actual body mass plus the added mass is called the 
virtual, induced, or total mass. Similarly, there are added moments of 
inertia terms which accompany the body’s angular acceleration. Conse- 
quently, for a long flat plate of dimensions such as 

- : : : :J 

2 a 

the added mass of the entrained water for motion perpendicular to the plane 

W 

of the cradle is defined as = pma^, or in our case a = ^ CLeptn = L so 
= ^PttW^L. The total mass is then actual mass plus added mass 

Mp = M( actual) + M (added) 

. Wp = 508 + JtipTrW^L. 

The added roll moment of inertia due to the entrained water is defined as 
or in our case 1^2 = l/SpTra^*, 

Ia2 = 1/8 ptt(|)‘*L. 

The total roll moment of inertia is then 

MVJ 2 1 W** 

^2 “ 12 8 16 



Similarly, for pitch it is 



^ ML2 ^ 1 „ 

I = + -- pTT — VJ. 

1 12 8 16 
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r * ^ ^ ^**^1 >■» »^ r *. 

These forces consisting of buoyancy and weight are constant and thus 
do not enter into the equation in any way. 

Overview 

The equations developed to this point are nonlinear due to the drag 
terms but fortunately can be solved by various numerical methods of inte- 
gration. The equations, arranged in state variable form are shown in 
Figure 9 • It is important to note that the drag terms shown are extremely 
complex and are functions of five of the six state variables. More specif- 
ically, the heave, pitch, and roll drag forces are functions of the follovr- 
ing: y(2), vertical velocity of cradle; yCS), the cradle pitch angle; y (^)5 
the cradle pitch angular velocity; y(5)> the cradle roll angle; and y(6), 

.w— - ^ ^ , Z -• '-'J--*- 

drags are different functions of these variables. 

The forms of these functions have been explained in previous ’par a- 



’ graphs. 
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-DRAG (heave) 

f(y(2) ,y(3) .y(^) ,y(5hy(6)) 
finite element 



(KiL+K 2L-K3L-KI*L 



2y^ 



(-KiW-K 2*'^-K3 W+Kj^W) 

2i*Lp 



0 



(y . y.j.y j.y “DRAG ( ? i t Ch ) 

(KjL-.i. .<1 -aP r{y(2),y(3),y(l,),y(5),y(6)) 



finite elcner.t 



(K3L2 +KuL^+K 3L?+?C2I 2) 






( ♦K3LW- K^LW+K^ LW-K2WL) 
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chapter III. 

DEVELGPMEIiT OF LIIIEAR APPROX I MAT I OHS 

As discussed previously, the equations as developed are nonlinear. 
Unfortunately, this fact results in the necessity of a brute force ap- 
proach to the solution. That is, the veil developed test and analysis 
methods vhich apply to systems of linear equations cannot be used to draw 
general conclusions about a nonlinear system. Consequently, linear ap- 
proximations for the drag terms would be helpful only if they were shown 
to give reasonable results when compared to the exact approach. 

Ideally, the nonlinear drag terms, which now are functions of five 
of the six state variables , should be reduced to functions of only one 
velocity term. Actually, it would be reasonable to expect the majority 
OX a. puruj.cuxa.r ux*a^ ueriii Lu ue udut5t=;vl Llit, 

with that term. That is, heave drag should mainly be caused by heave 
velocity, pitch drag by pitch angular velocity, and similarly, roll drag 
by roll angular velocity. 

Therefore, using this as an initial approach, the single variable 
heave drag equation might be: 

Drag (heave) = pCj^WLy^ 

The equation for pitch drag would be somewhat more complicated and would 
be calculated as shown in Figure 10. 
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Similarly the equation for roll drag is 
DRAG (roll) = 



pCpWL^^ 



b4 

As yet, even these equations are nonlinear due to their squared velocity terms. 

Several approximations to velocity squared damping have been used in the 
past. One of the commonly used electrical approximations consists of equat- 
ing the energy dissipated by a current squared resistor to the energy dissi- 
pated by a linear resistor when both are driven by a sinusoidal current 

source. VThen this is done the term is in effect related to a constant 

8 

multiplied hy a single power of I. The constant which results is — Iq“]_ 
where Iq is the raaxira\ua current value and is the driving frequency. 

Therefore, if we consider the sea as a velocity squared damper, a rea- 



sonable approximation to heave drag would be 
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where Eaxinujri velocity of y. If this approach is applied to 

pitch and roll the follov;ing equations result: 

l4.0pCpLWYouiy 

Drag (heave) - r 

3:t 

pCpLH/6oa)ie 

Drag (pitch) = pit it 

pCj5W‘*L<I>ot^2'^ 

Drag (roll) =' ~ 

To check the validity of this linearization approach as well as the ac- 
curacy of the 6h element method, a comparison vas made with values obtained 
from the exact drag equations shown in Figure 10. Since these exact equa- 
tions are only valid for cases where single velocity terms are used, the com- 
parison was done for a case where heave drag is only a function of heave 

drag is only a function of roll angular velocity. In other words, simple 
roll, pitch, and heave drag were tested independently. 

The results of the comparison shown in Figure 11 indicate clearly how 
the finite element method compares almost identically with the exact values 
as computed with the equations of Figure 10. The one difference, that being 
that the linear approximation gave lower peak values as they appeared more 
sinusoidal in shape was anticipated because we are equating energy dissipa- 
tion rather than peak values. 

Additionally, to assure that the one teim approximation will give the 
desired drag force at the proper time, a comparison was done between actual 
and linear drag values which would be created during normal cradle oscilla- 
tions. Figures 12 through ik compare drag values with Lulu exhibiting a 
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Figure 11 
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ten neecnd ?.0^ roll and pitch, and a h.2^ ft. heave amplitude. The compari- 
son chovs that for this particular driving frequency and amplitude the linear 
approximation is in phase with the actual values, yet the amplitudes do not 
correspond. Another comparison at twice the frequency of the first confirmed 
that while a single variable drag term approximates the timing of the drag 
force well, it does not predict the amplitude. 

In an attempt to better estimate the drag forces several two and three 
term approximations were investigated but the results showed that, in gener- 
al, complete confusion reigned. Even worse, the drag forces were no longer 
in phase with the actual values. 

The incorrect drag amplitudes result mainly from the terms in the coef- 
ficient of proportion of the linear drag term. In particular, for heave 






nn - the freauency of the cradle. 



^11 ^ 

and Yq- > the maximum displacement of the cradle. The cradle frequency 
can be assumed to be the same as the driving frequency, but Yq , the maximum 
cradle amplitude, cannot be assumed to be the same as the maximiam driving 
amplitude. In actuality, Yq of the cradle is the maximum displacement of the 
cradle as given by the equations of motion. However, that value of Yq de- 
pends on the chosen value for the drag term, which in turn is a function of 
Yq. Basically we have an iterative problem where values for the drag terms 
must be initially assumed and comparisons must be made with the actual non- 
linear responses. 

To provide a baseline for comparison a brute force analysis was con- 
ducted on the nonlinear equations to determine the actual response or gain 
of the system at various input frequencies. Ten separate frequencies were 
analysed and the steady state gains were plotted versus frequency. The same 
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anaiysii; -was reaeatoci but this tin.e using the singie term linear drag ap- 
proximations in place of the accepted 6b element drag values. This analysis 
was also conducted at the ten different driving frequencies with three dif- 
ferent driving amplitudes . 

Figures 15 and l6 show clearly the differences between the linear and 
nonlinear responses at various driving frequencies. The roll response curve 
is ommitted since it falls just short of being on top of the pitch curves. 
These figures show that the linear approximations peaJc at the same frequen- 
cies as do the actual response curves but the gains vary widely depending 
on the input amplitude. In general, the linearized response curves show a 
condition of lesser damping than is actually the case. However, as the driv- 
ing amplitudes go up, the apparent damping increases due to the presence of 
the larger value of driving amplitude (1^)) in the coefficient of proportion. 

- Since the main value of such an analysis is the prediction of the 

resonance peaks , and since the linear and nonlinear responses are almost 

w u 

identical excent in the rangeO.5 7T to 2.0 — , it is valuable to correct the 

n 

linear response curves by increasing the apparent damping until the linear 
responses con-form to the nonlinear responses. This can be done by multiply- 
ing the previously discussed linearized drag terms by a constant depending 
on the magnitude of the driving amplitude . 

Comparing the linear response curves to the simple single degree of 

freedom system shown in Figure ITj we find that the gain of such a system 

2'q C 

at resonance is given by the equation — vrhere — is the damping factor. 

Cc Cc 

Using this equation we find that the nonlinear heave resonance peak compares 



well to a linear system with a damping factor of 0.193. Similarly the pitch 
response compares well with a damping factor of 0.16T. After finding the 
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indicated linear values of daxnping factors for the linear response curves 
with various input amplitudes, a plot of drag multipliers versus amplitude 
input was made. This plot, shown in Figure l8, indicates the factor by 
which the drag terms must be multiplied in order to make the linear response 
correspond closely to the nonlinear response for each input amplitude. For 
example, for our normal h.23 ft. heave and 10° roll and pitch amplitude, the 
heave drag term would need to be multiplied by 2.J, the pitch drag term by 
and the roll drag term by U.O. 

Slimming up, linearization of the nonlinear drag terms is a complex mat- 
ter. As in all problem.s in the physical sciences it is never feasible to 
express such linearizations with absolute certainty and with the inclusion 
of all possible variables. In fact, to treat any physical problem in a 
mathematical way it is necessary to first make assumptiens about the system's 
behavior that may not be entirely correct. In this case drag is assumed to 
be linearly proportional to velocity, but, as is later discovered,, it is also 
related to the driving amplitude and therefore appropriate correction factors 
must be applied. In spite of the difficulty involved, the value of this 
linearization is to provide the coefficients for the linear A matrix which 



can then be operated on to find the general trends of the system. 
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Figure l8 
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CHAPTER IV. 

SOLUTION OF EQUATIONS 

The complete solution of differential equations may be split into two 
parts . a steady-state solution and a transient solution which dies out at 
a rate depending on the degree of damping, mass, and time. Our main interest 
here is with the steady-state solution. 

Time Domain Analysis 

The tim.e domain solution of the equations of motion was found by the 
use of numerical integration. The integration was done using the newly in- 
stalled Interdata computer facility at M.I.T.'s Mechanical Engineering De- 
partment. The same department's program DYSYS, (DYnamic SYstem Simulation), 
was used xui' the lIiuc uumaj-u analvois of tl'.c wCnlmccr ccuct— cnc . DYSY? 
solves a set of first order linear or nonlinear differential equations and 
can be particularly useful for obtaining transient as well as steady-state 
solutions to nonlinear time varying problems. DYSYS uses a fourth order 
Runge-Kutta integration technique to obtain an approximate solution to this 
equation — 

Given Y« , the value of the Y vector at some value of t, say t , then 
^ .... 0 

the approximate value of the y vector at tQ + At is 

At 

Y + At = Yq + -^Kq + 2 % + 2K2 + K3) , 
the K's being equal to 



Kq = r'CtQ, Yq) 

^1 = Yq + 35 Kq) 



:! 
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l<2 = f'<to + V() + HKj) 

K3 = F(to + At, Yq + K2) 

The Runge-Kutta !!ethod is derived from a Taylor Series Expansion about the 
point tQ and . Theoretically, the method will converge to an exact solu- 
tion as t approaches zero... The method has been shown to be stable, self- 
starting and to have a built-in error check. 

Since the program is set up for simultaneous first order differential 
equations, the three second-order equations must be expressed as six first- 
order equations. This can be done by letting 

y = y(i) 
y(2) = y(l) ■= y 
y(2) = y 

Similarly , 

yd) = y 
y(2) = y 
y(3) = 6 
y(4) = 9 
y(5) = (f) 
y(6) = ({), 

and hence y(l) through y(6) are the six state variables which completely 
define the system. 

The computer program v;hich was written to effect the solution is given 
in the appendix. The program has as subroutines EQSYM, which indicates 
"Equation Simulation’’ as well as three separate drag subroutines which must 
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be called each time step and calculates the drag at each point in the 
program. A time interval of 0.2 seconds v;as selected as a time step and 
gave reasonable results. 

Frequency Domain Analysis 

Since not only specific solutions for given parameters were desired, 
a frequency domain analysis resulting in Bode plots was \mdertaken. Bode 
plots clearly illustrate the relative stability of a system, with gain and 
phase margins often being defined in terms of Bode plots. 

This analysis was conducted \:sing the M.I.T. Mechanical Engineering 
Department computer program MBODE. This program calculates frequency do- 
main gain and phase plots for linear dynamic systems (with constant coef- 
ficients) expressed in matrix form such as 

X = Ax + Bu 

y = Cx + Du 

where 

A = the system A matrix 

B = the input matrix 

C = the output matrix (usually diagonal ones or 
appropriate weighting factors) 

D = the direct coupled output matrix (usually zero). 

The fact that MBODE only operates on linear constant coefficient equa- 
tions is the reason we were forced to develop linear approximations for the 
nonlinear drag terms discussed .in section III. The Bode plots which 



resulted from this method of analysis show clearly the properties of the 
system and these are presented in section V. 
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CHAPTER V. 



DISCUSSION OF RESULTS 



Time Domain Response 

For this analysis numerous trials were run at various damping levels, 
input frequencies, input amplitudes, and spring constants. The majority 

of the trials were run at the "baseline conditions of h.23 ft. ship heave, 

0 Ibf ^ 

10 roll and pitch, and 20 ~~ spring constant with a 10% deviation among 

X o 

the four springs. The 10% deviation* is the method by which the nonideal 
nature of the springs was modeled and is the vehicle for energy transfer 
between the different motion modes. 

In a discussion with J. Mavor of W.H.O.I. it was learned that the pre- 
viously mentioned amplitudes with a 5-6 second period could be considered 

■; c»4-n-4-c^ -CrMiv* o o o TTt IQ cVir\T.rc +‘hnc 1 i Q+.n o 

situation of heave, pitch and roll and indicates that the transmitted ship 

heave is reduced from It. 25 ft. to about 5 inches while the roll arid pitch 

are reduced from 10° to 2°. A total effective spring constant of 80-100 
Ibf 

is required to accomplish this decoupling; however, the analysis il- 
lustrates the fact that ship-cradle motion can be uncoupled if soft springs 
are used. 

In Figure 20 the pitch and roll amplitude have been doubled and the 
ship has a slow 10-12 second period to give a "worst case" example of pre- 
dicted cradle response. Although this decoupling is not as good as with 
the shorter period, the heave has still been reduced to less than one foot 
while the roll and pitch have been reduced from 20° to' 10° and T° respec- 



tively . 
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Frequency Donain Response 



Since specific response curves do not reveal much about the general 
nature of the system ^ a frequency domain analysis was conducted in order to 
get more descriptive information on how the system will react to various 
input frequencies. Figure 21 is a plot of the system’s gain response at 
baseline conditions of amplitude and spring constant but at ten different 
input frequencies. It is , in effect, a Bode plot of the nonlinear case. 

The plot shows clearly the three reasonances of the system, the leftmost 
peak being the heave resonance and the rightmost peaks being pitch and roll 
resonances, the roll resonance being the smaller of the two. 

This figure required ten separate solutions of the nonlinear -equa- 
tions of motion involving great expenditures of computer time. Although 
the plot shows the resnonse of the system under fixed conditions of spring 
constant, mass, and driving amplitude, it is clearly not a desirable way in 
which to find how the system responds when differing sea conditions exist 
or when the system parameters are modified. It is here that the value of 
the linearization process comes into full view. 

In Chapter III it was noted that the nonlinear drag terms could be 

linearized once the driving amplitude was known. Again using the U.25 ft. 

, o 

heave with 10 roll and pitch as the baseline we can find the linear drag 
approximation by referring to Figure l8. Chapter III, and obtaining the 
proper multiplier. Once these drag coefficients are set. Bode plots may 
be initiated which show the response of all outputs versus all inputs over 
a wide frequency range. Additionally, the effects of .variations in spring 
constant or mass may be studied with only minor modifications to the now 



linearized A matrix. 
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Figures 2.2 through 25 show the system's response across a range of 
.01 to 1.0 Hertz or a roll pitch and heave period of 1.0 to 100.0 seconds. 
HjC outputs for heave versus ship pitch and roll show that the ship's 
pitching and rolling couples energy into heave motion which then peaks at 
heave natural frequency as well as the pitch and roll natural frequency. 
This fact indicates the dangers of designing a decoupling system to de- 
couple only heave. It is clear from the plot that a decoupling system de- 
signed only for heave may well excite the pitch and roll modes at their 
natural frequencies. This fact was overlooked ty single degree of freedom 
analysis . 

Several other conclusions may he drawn from Figures 22 through 25. 
While a 20 Ibf/ft was chosen as a baseline for the spring constants and 



was Used for Figures 22 and 23, it would be desirable to predict easily 
what would happen with higher or lower spring constants. Figures 2h and 
25 show the response of the system with a spring constant of four -times 
that of Figures 22 and 23* It is not surprising to see that the resonant 
frequencies are shifted upward by a factor of approximately two. This in- 
dicates that in spite of the nonlinearity of the parent system, the system 
obeys the basic laws for calculating natural frequencies of translational 
and rotational systems. Being more specific, heave natural frequency is 



/Keff 

proportional to • , while pitch and roll natural frequencies are 

proportional to ^57 and respectively. 

II 12 



As a result of the above, we have information on the system's behavior 
over a wide range of frequencies without the need of calculating a Bode plot 
for each value of K. In effect, using K = 20 Ibf/ft as the starting point. 



we need only shift the response peaks up or down the frequency scale by 



This sane method of peak shifting can be applied to variations of mass or 
moment of inertia in a like manner, but in this case the shift will be 



Summing up, this analysis has shown that motions of a subm.erged cradle 
can, in fact, be nearly decoupled from surface ship motions if a soft ef- 
fective spring constant is provided at the four cradle corners. In addi- 
tion, the analysis has shown that the critical design consideration is the 
pitch and roll natural frequency rather than the natural frequency of 
heave, *= ° voQuit . careful thought must be given to the selection of 
spring constants since the equations of motion are coupled and are capable 
of transferring energy from one mode to another and hence exciting other 



a factor of the square root of K 




K20 



proportional to 
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modes at their natural frequencies. 
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CHAPTER VI. 

APPLICATION POSSIBILITIES 



It is not the major objective of this thesis to offer a detailed de- 
sign of a decoupling system. It is enlightening, however, to discuss 
several approaches that could be applied in the design of such a system. 
The primary consideration here is to obtain a soft spring constant and 
still support a large static load. 

From the previous analysis it is apparent that the soft spring con- 
stant must be at the four cradle comers. Further review indicates a 
spring constant between 20 and 50 Ib^/ft would seem to be a reasonable 
selection for these comer springs. 

At this noir.t the’^'^' njo rvoes of approaches to coiisider. The 
first and simplest approach is the passive system where we try to create 
a spring with the proper softness so that without any outside control the 
springs will decouple the ship-cradle motions. This approach can be 
likened to vibration isolation of a motor or a vibrating mass. It is sim 
pie and reasonably foolproof. 

The alternative method is of course a more complex active system 
where some sort of control is used to actively uncouple the two motions. 
This control can be in the form of position, velocity, acceleration, inte 
gral or differential feedback, horizontal sensing transducers or mechani- 
cal linkages. Consequently, unlike passive systems, these types of sys- 
tems can be timed to changing sea conditions or conditions of loading, 
thereby giving them a possible advantage. 
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E’assive Syntoms 

To create a passive system for the cradle v;ould require consideration 
of several methods of achieving the desired soft spring constant. A first 
choice for component selection would probably be simple metallic springs 
vjhich deform in accordance with Hooke's Law. Unfortunately, the static 
cradle weight of about 16 ,000 pounds would stretch the four 50 Ib^r/ft 
springs to about 80 feet before reaching equilibrium. Clearly, this would 
be unacceptable for if coil springs were used they would necessarily re- 
semble long "slinkys’’ and v/ould prove completely unmanageable. Leaf springs 
vjould similarly be of little use as they too would deflect excessively in 
the static mode. 

Torsional springs would possibly work if they were rigged through a 
complex geax-ing system v/hich woi.md up as load was applied; ho'wever, they 
too would require a large static deflection under the simple cradle load 
and would involve complex machinery and cable drums. Such a mechanical sys- 
tem would probably be more of a liability than an asset in an ocean environ- 
ment. 

Another possible alternative might be to take advantage of the non- 
linear stress strain properties of elastomers. Vulcsinized rubber will 
stretch to nearly ten times its own length and support up to ten tons per 
square inch. Unfortunately, the resistance to deformation usually increases 
in greater proportion than does the extension. If, by special treatment, an 
elastomer with the opposite properties were developed, it might be a logical 
choice to use as a spring. However, a short literature search turned up no 
suggestion that such a material exists. 



-59- 



A more roanonable method of producing such an effective spring constant 
would be with a hydraulic or pneumatic device. A device such as a hydraulic 
or pneumatic cylinder has the advantage that it can be biased to an initial 
static load and then oscillate about that equilibrium position with an ef- 
fective spring constant which depends on pressure , temperature , properties 
of the fluid, and dimensions of the system. 

To illustrate such a system consider Figure 26. 
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This system would have a spring constant which would depend heavily on k, 
the compressibility of the fluid. It is evident that a fluid with a high 
compressibility, such as air, would result in a softer overall spring con- 
stant which is desirable in this case. Additionally, a decoupling system 
with air springs would be capable of being set to an initial static force 
equal to that of the cradle weight. This might be accomplished with a sim- 
ple pressure regulating device. Following initial set, the cradle could 
then oscillate about equilibrium with a soft effective spring constant de- 
pending on the previously mentioned parameters. 
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Before considering this simple approach there are two important con- 
siderations to be aware of. First, the oscillating piston must look exact- 
ly like a soft spring to the disturbing ship motion. That is , the differ- 
ence in piston force from equilibrium to the extremes of piston travel must 
compare well with what would be seen if a real metallic spring were being 
used. In other words, for a ten foot pneumatic cylinder with a normally 
centered piston (capable of ± 5 foot excursions), the fully extended piston 
force must be equal to static cradle weight plus K (effective) * 5 feet. 

For example, if the desired spring constant were 50 Ibf/ft, then the piston 
could only exert 250 Ib^ more than the cradle weight at its maaimum stroke 
of five feet . Quickly calculating the change in pressure for full compres- 
sion using PqV 2 ^ = obvious that if the pneumatic system were 

dead ended, the reduction in volume to i/5 of its previous value would re- 
sult in a fivefold increase in pressure. 

Consequently, the first consideration in such a design is the fact that 
a plenum chamber would have to be used in conjunction with such a system in 
order to have the desired results. To illustrate the size of the required 
plenum chamber, we can again use = ?2^2 (if = T 2 )» Assuming a base- 

line piston area of 1 ft.^, a static cradle load of l»000 Ib^ per spring 
would require a piston pressure difference of 2T8 psi in order to maintain 
the cradle at equilibrium. Calculating the plenum size required with the 
following conditions 

AV at full stroke = 5 ft . 3 
AF at full stroke = 250 Ib^ 

PlVi = P 2 V 2 or 278 = (278 + 1.7) (V^^ - 5 ft^) 

= 815 ft^ 
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ve unfortur:ately find it is a chamber far exceeding the capacity of the 
cradle to easily carry. 

Hov/ever , further calculations show that the size of the required plenum 
will go down as working pressure goes up. For example, 
assume = 2000 psi 

Piston area = 2 in.^ 

or 2000 = (2000 + 125) (V^ - .069) 

= 9. It ft.^ 

and is a volume vhich is well within the carrying capacity of the cradle. 

The second major consideration in passive system design is the problem 
of centering the system. With a pne\imatic spring such as is sketched in 
Figure 27, the piston and plenum might be pressurized on the surface until 
the Tiiston was positioned at midstroke and supporting the chcuile's weight. 




Figure 27- 
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Tiiis initial setting, although correct on the surface, would be incorrect 
after the effects of hydrostatic pressure, temperature, and leakage had taken 
place. This fact would necessitate the constant bleeding in of make-up gas. 
In effect, some active control would be necessary to make this passive design 
workable . 

In addition to the aforementioned major design considerations, when us- 
ing a gas, with its higher compressibility, there are serious technical prob- 
lems to be overcome. These problems involve containment of a gas under high 
pressure, sealing, and danger of explosion, all of which may lead to bulky 
and expensive design. 

Throughout this entire first section only passive decoupling systems are 
discussed. These systems all have the same characteristics in that they os- 
rjinoti, ftbont an enuijibrium uositicn not deuermined actively by sensors end 
feedback, but by cradle weight, system dimensions, and initial conditions. 



Active Systems 

Systems which employ active feedback for control are called active sys- 
tems. Typical examples of these systems include automobile power steering, 
hydraulic positioning devices, and automated milling machine tables. These 
active systems are all characterized by a mechanical, electrical, or hy- 
draulic system by which information on the output is fed back to control the 
actions of the input. This feedback control can be designed to provide fast 
or slow response, maximum or minimum overshoot or undershoot, and nearly any 
natural frequency within the constraints imposed by the parameters of the 
system. 

The simplest form of an electrohydraulic position control device is 
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shovm below in Figure 28. 




Figure 28 . 

The feeUbacK loop consists oi sensors whica leeu oac^ posioioa, veloc- 
ity, and acceleration information which improve the system's response to 
disturbances at the input. Position feedback often is sufficient for system 
control but does not provide the tight feedback loop that results when 
velocity and acceleration are also used. 

A typical response curve for such a system is shown in Figure 29. 




Figure 29- 






That is, the system voulcl tend to rise to the desired level, overshoot it, 
and then oscillo.te about that level with a decaying oscillation depending 
on the parameters of the system. It is this oscillation which might well be 
used to create a virtual soft spring constant at the constant load. 

A Proportional Control Approach 

As previously stated, feeding back all three variables v;ould tend to 
make the system have less overshoot and oscillation. Since we are attempt- 
ing to utilize the system's natural oscillations, we would be more success- 
ful using only a position feedback system which, while not enhancing oscil- 
lations, does not discourage them strongly. 

To determine the feasibility of such a system and to get some design 
parameters for guidance, the system equations are determined below. 

First, looking at the spool valve and assuming the hydraulic fluid is 
water , 

Valve flow — 




where = flow to load 

Kqx = flow sensitivity or valve gain 



Kqp = negative slope of valve characteristic curve at 



constant value of x. See Figure 30. 




Figure 30 . 
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For an ideal valve the slope, is very small and can be assumed as 

Such a valve is called a flow control valve. 

Valve forces — 

-KrX + F. = 0 

where Kp = valve control spring 
= valve force 

Next, considering the piston, 

Piston flow — 

= # * P„ !s (ICe 

where : 

C 2 = piston leakage 

Kg = compliance of piston ' " 

V = volume of piston chamber 

= bulk modulus of fluid 
m 

Piston forces — 

M" = P A + V^^. 
y m t 

A = area of piston 
M = mass of load 
.W = weight of static load 

Combining all four equations , 

m 

and : ^^y = P A + W^. . 

in ^ 

Finally, 

^ 0'l(Ke I )) ^ 



CpMy 

^ Ay = Kg^X 



zero . , 
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and usins LaPlace trancforms , 
y(s) = 



Ko,.X(s) 



(s2-M.(Ke + + s + A) 

- - ^ Pm A 



2A 



Such a system can be described in terms of poles and zeros and valuable 
information can be obtained with an S plane plot of the position of these 
poles and zeros. Our case results in three poles, one at the origin and a 
pair of complex conjugate poles. The conjugate pole pair indicates the 
system is oscillatory while the pole at the origin indicates that the system 
may not reach the desired new setpoint or that it will have a steady state 
error. However, this is of little concern to us here. Plotting the zeros 



in Figure 31, 




Figure 31. 



The natural frequency of interest is defined by the length of the 
phazor connecting the origin to a complex pole. Basically, if we put the 
transfer function in the form 



T(s) 



const 

s(s + + j3j^)(s + - jB^) 



we can adjust the natural frequency of the system by adjusting the coef- 
ficients which control the position of the complex conjugate poles. 






I 





In our case, 
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calling = ^2 



Mp 



= An 



and rewriting 



y(s) = 



= A, 



const 



s{A^s^ + Aq^s + Aq) 



The natural freo[uency of the system is defined by the roots of the 



equation 



or the magnitude of 



or 



U An 



At _ / ,A-\ ,2 . , . , , An , 

A2 A2 A2 



/w - 






If we assume leakage is equal to zero then ^2 = 0 and all terms with A^ 
go to zero, and the expression reduces to 

iHa 



0) = 
n 









MCK^ + 1) 
e e 

To get a natural frequency of 0.03HZ or about 0.2 the argument of the 
square root vould have to be about O. 5 O. 

Assuming a baseline piston size of 1 ft^ and V = 10 ft^ using a mass of 
2000 slugs and the bulk modulus of water 6= 320,000 psi, the equation gives 



- « • 








1 



0.50 



sec 



2 ft2 
2000 slugs (K„ + 



10 “ 



or, with proper units 



0.50 



32 ,000 



(Il4l4)(ll4l4)(2) 



) , 



sec2 " 2000 . lT28£inL) 



ft 



12in 320,000psi 



- 1497.00 



psi 



in^ 



Hence a conpliance of about 500 must be built into the system. Assviming 

an air accumulator with a working pressure of 278 psi and change of one psi 

PlVi = P 2 V 2 

278v^ = 279 (V^ - 500) 

= 80 ft^ , 

rmnr^'h ocrnTn q v<=‘T'V 1 pnn hllTkv. 

Fortunately, as was shown in the passive systems, this large accumulator 
may be made smaller by use of a higher working pressure. For example, 
assume P^ = 2000 psi 



Area = 2 in'” 



0.50 = 



2(1*) in'* 



2000 slugs 1ft / 20 

' *?on non ) 



12 in 



"e 320,000 



K = 0.096 ^ 

e psi 

or again, using an air accumulator at 2000 psi, one psi should change the 

3 



0.096 



in 

psi 



volume 
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2000 = 2001 (V^ - 0.096) 

192 in^ = 





which is of reasonable size and is easy to build. 

Actually, a working pressure somewhere beWeen the limits discussed 
here w’ould be most appropriate in order to reduce piston size and weight 
as well as reduce the plenxim size. In short, there is probably an optimum 
selection of working pressure, piston area, and hydraulic fluid. A no- 
tional system employing this approach is shown below in Figure 32. 







Figure 32. 
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A Proportional-Integral Control Approach 



By increasing the level of control is include integral as veil as 
proportional control we may be able to improve the system and possibly re- 
duce the size of the plenum chamber. 

Modeling the system as shown below the following equations result. 



Figure 33. 




Combining equations (l), (2), and (3) 



My = - - Fowsinoit - Dy 



but 



p = _i_ + ^ 

^m Ket AKet ’ 



Area = A 



Pjjj = Load pressure 

Ket “ Total fluid plus 

cylinder compliance 

0^ = Flow to load 
FoSinoit = Disturbance 

R = Proportional con- 
trol consiam. 

S = Integral control 
constant 

D = Damping constant 

. • ^ 

(1) y - Kg^Pjjj - ^ 

(2) My = -Pjj^A-FQsimot-Dy 

(3) Qm = Ry+s/ydt-c(t) . 



■ Plugging in 



my = 



- A(j^ + sfe') - Fowcosiot - D^^ 




I 

i 

I 
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also , 

Qjj, = Ry + Sy - C(t) 



or finally , 



^ " K,t 



R • 

y + — = C( t )+Fow^sino)t . 



K 



et 



Using LaPlace transforr.s, 



y(s) = 



C(t) + Fow^sino)t 



(MS'* + Ds3 + -^ S2 + 5§. + _S 



K, 



et 



K, 



et ^et 



V/here C(t) is the input setpoint (time rate of change) and Fov^sinmt is the 
disturbance 5 the characteristic equation is then 



MS'* + DS^ 





S 

Ket 



0 



or 






■no j n 

+ ;r7‘ 

M KetM 



■’"‘C 

ic^ 



i^tM 



u . 



This equation can then be analyzed with standard procedures and proper co- 
efficients for R, S, and Kg-j. can be found. 

Some general comments can be made about proportional-integral versus 
simple proportional control. Addition of integral control allows the sys- 
tem to respond to a change in set point with zero steady state error while 
it does not change the basic natural frequency of the system. Moreover, 
this control element provides an additional independent constant which can 
be used to more accurately tune the system for maximvun decoupling effect. 
The availability of this additional constant may well offset the complexity 
at^ded with an integral feedback unit. 

One fact to be aware of with the P.I. control discussed above is that 
the cradle has been shown to exhibit less than critical damping. Hence, 
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the characterii:tic equation would probably become unstable upon the ap- 
plication of a reasonable gain and additional damping may be needed. This 
damping could be provided either by mechanical means or by additional feed- 
back such as velocity and acceleration. 

The determination of the optimal design for such a feedback control 
system is left for future work. 

Other Systems 

Basically, all other systems are combinations of the passive and active 
systems. Some system modifications which may be desirable include: 

(1) Inclusion of horizontal sensing valves which tend to channel fluid 
to the low sides of the cradle thus giving horizontal feedback to an active 
system. 

(2) Substitution of rotary air or fluid motors with gear racks to 
provide the vertical position control. 

( 3 ) The utilization of an active-passive system which employs a posi- 
tion control to center a passive piston or rack. 
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CHAPTER VII. 

COIJCLUSICHG AIID R.ECC:T lEHDATIOIlS 



As a result of the investigations carried out in this thesis, several 
conclusions have been arrived at. 

1. A four point passive decoupling scheme can substantially decrease 
ship heave, pitch, and roll trsinsmission to a submerged cradle. 

2. The effective spring constants required for decoupling fall in the 
range of 20-50 Ib^/ft. 

3. Each mode of oscillation will couple energy into the other two 
modes , hence the design of a decoupling system must take each natural fre- 
quency into consideration. 

Cradle hydrodynamic drag is a complex function of frequency, ampli- 
tude, and velocity and mus'*. De cnosen uo suj.o uae xiiLLj. V a.uuo._L o -L j-wi* 
is to be modeled. 

5. Conventional springs such as coil or leaf springs which will have a 
soft spring constant and also a large initial static load do not exist. 

6. Passive pneumatic springs, while meeting desired requirements, would 
require large plenum chambers in order to operate at safe and reasonable 
pressure levels. Additionally, these passive devices would be difficult to 
center in the cylinder. 

7. A logical approach for creation of soft springs with initial static 
loads would be an active feedback control system. A system, such as one 
utilizing proportional and integral feedback, could substantially decouple 
ship motions from the cradle, yet maintain the resultant piston oscillations 
about the center of the cylinder or set point. 
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In summary, both the theoretical and practical aspects 
of a cradle mounted decoupling system have been investigated. 
The value of this analysis is that it brings us to the point 
where we must decide whether the possible improvements in op- 
erating conditions warrant the considerable time and expense 
involved in making such a system come true. If the analysis 
had shown that a simple passive system v/ould work and was 
readily available, further tests and construction would be 
easy to justify. On the other hand, since an active system 
seems to be necessary, extreme caution should be exercised 
prior to construction to assure that; first, proper feedback 
coefficients are found; and second, the necessary control com- 
ponents are available and are readily adapraoie ro the oceeu 
environment. In addition to the above, care must be taken to 
insure that any resultant detailed design satisfys all safe- 
ty requirements which may be imposed by WHOI, the Navy, or 
any other regulatory bodys . 

Other means of decoupling ship-cradle motion have been, 
and are being evaluated. The author hopes that this work will 
enable WHOI to more clearly evaluate the cradle mounted system 
and how it compares to other methods being considered. 
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